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1. INTRODUCTION 

The control of the induction motor (IM) that has no speed sensor (sensorless) has the advantages 
such as low cost, high reliability, and saving space, reduces hardware, works well in chemical, dusty, 
decreases maintenance requirements. Due to the above advantages, there are now a number of methods for 
estimating the speed of induction motors without speed sensors studied such as the methods using the 
machine model is an estimative type of the open-loop [1], model reference adaptive system (MRAS) [2]-[7], 
the Luenberger observer [8]-[11], the sliding mode observer [12], the Kalman filter [13]-[18] (the first 
estimative types), and the methods which use machine models that are the estimator types using algorithms 
intelligent such as neural-networks [19]-[22], fuzzy-logic based on a control without the speedy sensor [23]-[25], 
(the second estimative types). Each method has its characteristics, advantages, and disadvantages, which will 
be presented below. 

The MRAS is based on comparing the two outputs of the two models: the first model (the reference 
model) contains no the rotor speed, the second model (the adaptive model) uses the speed to estimate the flux 
of the induction motor. The outputs of the two models are compared together to get the error value. The error 
is the input of an appropriate adaptation mechanism to generate the estimated rate fed back to the adaptive 
model. The benefit of this method is simplicity, fast processing ability but low accuracy. 
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For the reasons mentioned above, this paper presents the method of the speed estimation by 
extended Kalman filter instead. It belongs to both groups mentioned above, because it uses the Kalman filter 
optimized by the genetic algorithm (GA), which has more advantages such as good noise filtering and higher 
accuracy. Since the equation of state of the induction motor is nonlinear, we will perform the discretization, 
and thus in each small cycle that is considered linear to apply the Kalman algorithm. This is a recursive 
problem combined with adaptation to estimate the parameters in the equation of the state of an induction 
motor. A notable problem in the Kalman filter is determining the components in the covariance matrix Q, R. 
Traditionally, we would define these matrices by trial and error method, but this method is time-consuming, 
and its accuracy is not high. Therefore, this paper presents how to determine these matrices by GA, which 
belongs to the second group, that of intelligent algorithms. 

Except for the introduction and conclusion, this article includes four main sections. Section 1 is the 
modeling DFOC of the IM control without speed sensor, sections 2 is the speedy estimate method reference 
frame - model reference adaptive system (RF-MRAS), sections 3 is the speedy estimate method current 
based - model reference adaptive system (CB-MRAS), sections 4 includes details to implement the speedy 
estimate of the EKF: design of proposed EKF algorithm, determination of components of matrix Q, R by trial and 
error method, and GA algorithm. Finally, the simulation results of three methods in two cases: the resistance of 
rotor is constant and change during operation process. 

The simulation results showed that the speed response of the Kalman filter optimized by GA is 
better than the RF-MRAS, CB-MRAS methods under different working conditions, when the rotor resistance 
varies with temperature. This also means that the GA-Kalman filter estimation method is less dependent on the 
system parameters than other methods. The studies covered in this introduction are presented in the next section. 


2. THE MODELING DFOC OF INDUCTION MOTOR 

From the equation system of an induction motor, the sensorless speed direct field-oriented control 
(SS-DFOC) for induction motor drive is constructed as following Figure | [26]-[29]. In this model, the proposed 
genetic algorithm - extended Kalman filter (GA-EKF) speedy estimate block with input parameters of voltage and 
current collected from the system through voltage and current sensors. The outputs are the estimated speed and 
flux current. 


Figure 1. The structure of the SS-DFOC method 


The vector control (FOC) includes controller of the stator currents characterized by a vector. This 
control is based on the design that converts a three-phase time and speed system into a two co-ordinate (x and y 
co-ordinate) time invariable system. This design lead to a configuration similar to a DC machine control. Vector 
controlled machines need two constants as input references: the torque cureent component represented by rotor 
speed wọ and the flux current component represented by Im [26], [27]. 


3. ESTIMATE SPEED OF INDUCTION MOTOR USING MODEL RF-MRAS 
3.1. The mathematical equations of the RF-MRAS model 
The mathematical equations of the RF-MRAS model are as following [3]-[5]. 


Yra = =U (usa ~ Rsisa)dt > dLsisa) > Wrp = UJ (use = Rgisg Jat 7 dLsisp) (1) 
Dra = (2 İsa — 7 Pra = Wore) dt; Drp = Ie isg — 7 rp + WoPra) dt (2) 
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This model only cares about the motor speed, so the value of the stator resistor is known in the simulation. 
The error signal is represented below by the expression (with Kp > 0, K; > 0). 


$= DraWrp = Deere: @o = Kpé + Kı IN édt (3) 


Where the state variables isg, isp, Usa, UspWra Weg, Wo are stator currents, stator voltages, rotor fluxes, and 
the rotor speed of IM in [a — f] system, the parameters with a hat on the top are the estimated values. 


3.2. The RF-MRAS model 

The RF-MRAS for the speed estimation has diagram such as Figure 2 [6]. The model consists of two 
blocks: the adaptive block, and the reference block. Two these blocks are compared with each other, the difference 
in values of two these blocks is fed back to the adaptive block to correct the desired signal for accuracy. 
The current and voltage signals are specified in the model. 


Angular 


velocity 
Adaptive 
Mechanism 


Figure 2. The diagram of speed estimation using the RF-MRAS model 


4. ESTIMATE SPEED OF INDUCTION MOTOR USING MODEL CB-MRAS 
4.1. The mathematical equations of the CB-MRAS model 
The Mathematical equations of the CB-MRAS Model is described (4) [6], [7]. 
z 1 A a z 3 1 * AF X 
lsa = 7 (Kiusa + Kora + Kz@rYRep T îs« Jdt; lsg = z S (Kiuse + K2Wrg — KzðRr Yra — îsg)dt (4) 


The rotor speed is calculated from the following expression with Kp > 0, K; > 0. 


; RAT ; a AROA t 
= Cisa ~ isa) Prg ~ (isg = isp) Pra» @, = Kpf + K; is dt (5) 
_ Lr : = Lm i p 1. es LsLy—L?, ` — Lrĝs Lm 
TES Calm? 2 ~ GA) K3 ~ Gy n= Cilm ’ as Lm t Tr (6) 


4.2. The CB-MRAS model 

From the above expressions, the CB-MRAS model is built below [6]. Figure 3 is the CB-MRAS model 
which is similar to the RF-MRAS model, but the block functions differ from the above model. The model has 
three input blocks: current model, voltage model, and current estimation block. The output has two blocks: 
the rotor speed adaptive block, and the stator resistance adaptive block. Error signals: rotor speed and stator 
resistance are fed back to correct the rotor speed and stator resistance values. 
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sar sp) Mo 


Figure 3. The diagram of speed estimation using CB-MRAS model 
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5. ESTIMATE THE SPEED OF IM USING EXTENDED KALMAN FILTER OPTIMIZED BY 
GA ALGORITHM 
In this section, the first presents the discrete form model of the equations of state. The second is 
the algorithm to the EKF, how to determine the expressions in the algorithm, the Q matrix in the system 
block, R in the measurement block, and limitations in determining the matrix Q, R by trial and error. Finally, 
the GA algorithm (a smart algorithm) is used to optimally determine the components in the Q, R matrix for 
estimating the speed of the IM in the most accurate method. 


5.1. Using EKF to estimate the speed of IM 

From the nonlinear equations with five state variables: isa, isg, Pra» Yre» Wo [26], [27], to be able to use 
the EKF recursively, They need to be transformed from the continuous state equations of the motor into a discrete 
form, and for each small cycle, these equations are considered linear. Then they are added into the system noise 
W and the measured noise V, so the EKF algorithm shown below [14]-[16] can be applied. 


Xn+1 = Án: Xn + Bn Un + Wa 
Yn+1 = Cn ' Xn+1 t Va (7) 


Here, the coefficient matrix is calculated according to the expression: 


An =e^ x I+ AT; B, = f) e“7Bdt ~ BT; C, =C 


T T T 
.(n+1) pry ary pee) wg? Yny = ag | Un = | (n) yp” 


Xn+1 = i SB lsa lsg = |Vsa SB 
1- Kr T 0 LaRR PLpwo 
Kı LRKI 2LRKi 
0 1-7 PLAWO LaRR IT o 1 07 
Kı 2LRK1 LRKI Kı 0 1 
An=| 2 1 Pog :B,=|9 <kce,=|o 0 8 
n= ee 0 Rees Dae, 0j; Dn = Ki) Cn = (8) 
: Ln pot 1 0 0 ay 
0 Te see 1- zT 0 0 0 0 0 
_ Pliny Plh an 
IE ARB aT ARa 0 0 1 
The covariance matrices Q and R of noises have the following form: 
Q, withn = L, R, withn = l 


Qn = Cov(w) = E[ww'] = { ; Ry = Cov(v) = E[vv'] = { 


0 otherwise 0 otherwise 


The matrices Q, R show that they have the form of diagonal matrices. The EKF algorithm is used for the 
speedy estimation model of the induction motors. 


Xn41 = AXn + Buy = f (xp, u”) (9) 
The linearization of nonlinear equation [26], [27] inplemented around the estimated value £; is as following 


— faxi u”) _ O(A&n-1+BUn-1) 


Fa Ox; Ox; (10) 
Pasi = R ÊF +Q (11) 
A(x") = aixi (12) 
Ih _ A(Cn(x{ xP 
H, = aan = a (13) 
Kn+1 = Pri CaP aah +R)* (14) 
Xn41 = Xn+1 + Kn410ns1 z CnXn+1) (15) 
Prat Si Prat = Kn+1CPn+1 (16) 
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Where the P is covariance matrix. The K is the Kalman filter gain. The x are state variables 
(parameters) of IM. The noise matrix Q, R, and the parameters with a hat are the estimated values. 
The beginning value of the covariance error matrix has the form below. 


4 0 0 0 0 10000 
0 a4, 0 0 0 0 01000 
Q=|0 0 à 0 OfR=|o ke=jo 0 1 0 0 
0 0 0 4 0 = 00010 
0 0 0 0 ae 00001 


To implement the EKF for estimating the rotor speed, current, and flux, the parameters Aj, 4; are 
appropriately chosen in the Q, R matrix. Traditionally, this is done by trial and error, which takes a long time but 
doesn’t produce very good results. For this reason, we will find the optimal solution for determining the parameters 
of the fit of Q and R by genetic algorithm. The algorithm will be presented in the following section. 


5.2. Using genetic algorithm to estimate parameters of Q, R of EKF 

To find the optimal parameters A;, 4;, we use a soft computational algorithm that is the GA [23], [30]. 
When the appropriate parameters are available, the motor speed estimation will be more accurate. The algorithm 
flowchart for finding parameters /;, 1; of the matrix Q, R is below. 

Figure 4 is the GA algorithm applied to this paper. In the first step of the algorithm, we randomly 
initialize the population (parameters Q, R) and encode the computation in the form of real numbers. This set 
of parameters will be transferred to the Matlab model to perform the simulation. Its results will be evaluated 
by the cost function, then perform genetic operations such as selection, crossover, and mutation. Thus, we will 
have a new set of values (usually better than the previous one) and continue to feed the model for testing. 
This process continues until the required number of generations reached. The result of the last run time is the best. 


Initial random population and encode to real number 


Generation=0 


Decode to 
The real number 


Evaluate fitness value for each 


chromosome 


Fitness value — — 


Gen > max 


generation 


Yes 


Optimal 
solution 


Figure 4. A flowchart of the genetic algorithm for the GA-EKF model 


6. SIMULATION RESULTS AND ANALYSIS 

The basic simulative parameters used for induction motors are as: P = 1.5 kW, Upc = 270 V, Pp =2, 
Rs = 2.1 Q, R, = 2.51 Q, Lm = 0.129 H, Ls = 0.137 H, L, = 0.137 H, J = 0.043 kg.m?. The three different 
speed levels simulated in this section are 100 rpm, 40 rpm, and —40 rpm. The simulation results in many 
cases are presented in the following sections. 
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6.1. Estimate speed of IM based on the model, RF-MRAS 

In the Figure 5 (the first figure of the simulation), the responses of the reference speed and the actual 
speed are shown. The solid line is the reference speed. The dashed line is the response of the actual speed. 
The actual speed follows the desired speed, but there is an overshoot. Figure 6 the comparison of Figure 6(a) 
the RF-MRAS estimated speed, the actual speed, and the reference speed and the relationship of Figure 6(b) 
the error between the estimated speed and the actual speed, we see that the speed error between the RF-MRAS 
estimated speed and the actual speed is large. Its peak value is 13 rpm, so this method is used in drives less 
than at present. 


— Reference Speed | 
120 ---Actual Rotor Speed 


-60 


-80 z | 
0 0.1 0.2 0.3 04 0.5 0.6 0.7 0.8 09 1 
Timeļfs] 


Figure 5. Reference and actual rotor speed of the IM drive 
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Figure 6. The comparison of (a) the RF-MRAS estimated speed and (b) the error between RF-MRAS 
estimated speed and actual speed 


6.2. Estimate speed of IM based on CB-MRAS model 

This section is the simulated results of the CB-RAS model. Figure 7 the evaluation of Figure 7(a) 
the CB-RAS estimated speed, the reference speed, and the actual speed and the analysis of Figure 7(b) 
the speedy error between the CB-RAS estimated model and the actual model, the deviation of the estimated 
speed of the CB-MRAS model is quite good. Its peak value is 2.7 rpm. 


Reference Speed 
120 Actual Rotor Speed 
CB-MRAS Estimation Speed 


Speed{rpm; 
Speed Difference 


(a) (b) 


Figure 7. The evaluation of (a) the CB-MRAS estimated speed and (b) the error between the CB-MRAS 
estimated speed and the actual speed 
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6.3. Estimate speed of IM based on EKF with defining matrices Q, R by trial and error 

First, the values A;, u; in the matrix Q, R which is chosen by trial and error are: A, = 9.23E-14; 
Az = 8.33E-14; A3 = 4.88E-14; A, = 6.85E-08; As = 1.46E-08; u = 9.13E-05; u, = 4.99E-06. With these 
values, the speedy response is shown in the Figure 8. The speed according to EKF algorithm is the dashed line. 
These values cause a divergence with a spiked error. This method takes a lot of time, but the results obtained are 
not optimal. We continue to probe by trial and error. The selected values are: A, = 8.74E-14; Az = 4.26E-14; 
Az = 1.69E-14; A, = 6.80E-14; A; = 3.26E-08; pW, = 1.79E-05; u, = 2.49E-05. With the new values for EKF 
method by trial and error, the speed response is displayed in the Figure 9. Its error is a small value that we can 
see in this figure with the dashed line. 


—Reference Speed 
— Actual Rotor Speed 
---EKF estimation Speed 


Figure 8. The reference, actual rotor, and EKF Figure 9. Reference, actual rotor, and new EKF 
estimated speed estimated speed 


6.4. Estimate speed of IM based on EKF with defining matrices Q, R by GA algorithm 

The two matrices Q and R are diagonal matrices of seven values that need to be determined. If the values 
(Ai A2, 43,44, As, H1, H2) are chosen by trial and error, it takes a long time, but the results described above are 
not very good. Because of this reason, the components in the matrix are found by the genetic algorithm to save 
time and get the best speedy estimation response. The Table | includes parameters of GA algorithm that we will 
use to find optimal values for two matrices Q and R. The range of values for variables is [le-18-0.1]. The cost 


Poe : ; 1 
function is calculated according to the error MSE is E = ~ di (Sreal_speed — Sestiiated-speea) - 


After performing the GA algorithm according to the diagram of Figure 4, we below collect results in 
the Table 2. Ten rows of the table correspond to 10 implementations of the GA algorithm. Columns 
Ài A2, A3, Aq, Às are of the Q matrix and columns p4, [2 are of the R matrix. The MSE column is the speed 
error corresponding to 10 runs of the algorithm. 

The best values after ten runs substituted into the matrix Q, R have the results shown below. Figure 10 
the analysis of Figure 10(a) the actual rotor speed, the GA_EKF speed, and the reference speed and the comparison 
of Figure 10(b) speed error between the GA_EKF model and the actual model. This speed error is the smallest 
compared to the two methods investigated above, the RF-MRAS and the CB-MRAS. Its peak value is less 
than 2 rpm, so this method can be used for IM drives. 

Table 3 is the results of the speedy error comparison. They consist of the RF-MRAS, the CB-MRAS, 
and the GA-EKF methods by the least squares method. The RF-MRAS Method has an error of 4.5502 rpm, 
the CB-MRAS method has an error of 0.1896 rpm, and the GA-EKF Method has an error of 0.0839. 


Table 1. The parameters of GA algorithm 
Parameters Values 
Population size 60 
The crossover probability 0.5 
The mutative probability 0.02 
Number of generations 10 


Simulation results of three methods, including graphs and errors evaluated according to MSE-standard, 
are presented above. In the proposed new method, the extended Kalman filter in which the parameters are 
optimized by the GA algorithm shows its advantages. These can conclude that the GA-EKF method is the best, 
although the CB-MRAS method was considered the best for a long time before. 
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Table 2. Convergence results after ten repetitions 


Index A, Az A3 Na As by H2 MSE 
1 4.82E-03 2.90E-10 2.90E-10 2.52E-09 5.84E-17 9.83E-03 5.84E-05 7.1586 
2 7.30E-15 1.08E-18 5.83E-08 5.84E-17 7.30E-15 0.042825 9.83E-14 5.0869 
3 2.98E-21 2.90E-10 4.28E-02 8.65E-01 4.87E-09 1.01E-04 9.65E-04 3.2586 
4 4.61E-02 7.82E-01 9.89E-01 5.31E-O1 7.50E-O1 1.12E-05 9.83E-04 2.1869 
5 4.84E-08 7.82E-07 9.89E-09 5.31E-07 7.50E-11  1.12E-04 2.92E-05 1.0821 
6 4.84E-08 7.82E-07 9.89E-09 5.31E-O7 7.50E-11  1.12E-04 2.92E-05 1.0821 
7 6.97E-11 2.90E-10 9.89E-05 8.65E-01 9.90E-01 1.01E-07  2.28E-06 0.5735 
8 2.26E-10 3.44E-12  2.26E-06 1.21E-12 9.83E-14 6.17E-06 6.17E-07 0.1286 
9 7.30E-15 3.85E-07 7.30E-15 7.61E-18 4.82E-10 3.85E-07 2.26E-06 0.0870 
10 2.65E-12 3.44E-16 2.65E-12 8.18E-18 2.26E-12 1.21E-07 1.21E-07 _ 0.0839 
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Figure 10. The analysis of (a) the EKF estimated speed and (b) the error between real speed and the EKF 
estimated speed 


Table 3. The error of the three methods 


n 
1 
= 2 
The methods E = — ) reat speed ET estimated_speed) 
1 


The evaluated results 


RF-MRAS 
CB-MRAS 
GA-EKF 


4.5502 
0.1896 
0.0839 


Good 


Very good 


The best 


6.5. The responses of the estimated speed in the cases of variable rotor resistance 


6.5.1. The resistance (R,) changes from 2.51 © to 2.0 Q 


During operation, the resistance of the rotor (R,) can be changed by temperature, this also means that 
the rotor time constant (T,) changes, it is a significant parameter in induction motors which affect the process of 
motor speed estimation. Therefore, its influence is considered in the three-speedy estimation methods mentioned 
above. In this section, the variable rotor resistance value and its influence on the speed estimation methods are 
investigated. The results of the three models are displayed in Figure 11. The variation of Figure 11(a) the rotor 
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resistance value from 2.51 Q to 2.0 Q and the comparison of Figure 11(b) the estimated speed of the RF-MRAS 
model, Figure 11(c) the estimated speed of the CB-MRAS model, and Figure 11(d) the estimated speed of 
the GA-EKF model, this is one of the cases in which the rotor resistance can change during operation. 

The simulated characteristics and results of all three methods which are mentioned above are 
summarized in Table 4. The first column of this table is the methods. The second column is the speed error. 
Finally, the third column is the evaluation result for each method. 

The results are obtained from the graphs and the calculation of the error according to the MSE 
standard, when the rotor resistance reduced by 20%. In third column, we see that the error of the GA-EKF 
method is the smallest (11.1203 rpm) and the speedy error of the RF-MRAS method is the largest. These 
show that the CB-MRAS method is better than RF-MRAS, but the EKF method is the best. 


Table 4. The error of the three methods (the Rr changes from 2.51Q to 2.09) 


1 
The methods £E = — ) reat speed — estimated speca)” The evaluated results 
1 


RF-MRAS 22.0262 Quite good 
CB-MRAS 12.4048 Good 
GA-EKF 11.1203 Good 


6.5.2. The resistance (R,.) changes from 2.51 Q to 3.0 Q 

Other changes in the rotor resistance are studied continually in this section. Now we will examine 
the case of increased resistance. From the simulation results, we obtain the graphs as well as the calculation 
of the speedy errors that will be presented in the following section. This section is similar to the one above. 
The results of all models are shown below. Figure 12 the change, Figure 12(a) of the rotor resistance value 
from 2.51 Q to 3.0 Q and the relationship of Figure 12(b) the estimated speed of the RF-MRAS model, 
Figure 12(c) the estimated speed of the CB-MRAS model, and Figure 12(d) the estimated speed of the GA-EKF 
model, this part compares the estimated speedy responses of all models in the change of the rotor resistance. 


— Reference Speed 
120) > Actual Rot peed | 
A ---RF-MRAS Estimation Speed 


— Reference Speed 


| 
zZ 


Actual Rotor Speed 120 —A peed 
N CB-MRAS Estimation Speed ---EKF estimation Speed 


Speed{rpm] 


(c) (d) 


Figure 11. The variation of (a) the change value of the rotor resistance, (b) the RF-MRAS model, 
(c) the CB-MRAS model, and (d) the GA-EKF model 
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All simulation results are listed in Table 5. Similar to Table 4, this Table 5 also summarizes the 
estimated results of all three methods: the RF-MRAS, the CB-MRAS, and the GA-EKF in the case of 
increased rotor resistance. The first column of this Table 5 the methods. The second column is the speed 
error. Finally, the third column is the evaluation results of each method. 

In this case, the rotor resistance increased to 20%. We evaluate the speed error according to the 
MSE standard. The RF-MRAS method has the largest error (13.5012 rpm), while the other two methods are 
quite good and approximately equal. These obtained results show that the GA-EKF method is equivalent to 
the CB-MRAS, and both these methods are better than the RF-MRAS method. 


Table 5. The error of the three methods (the R, changes from 2.51 Q to 3.0 Q) 


1 Th luated 
2 2 e evaluate 
The methods E= =) Caia = estimated_speed) 
1 


results 
RF-MRAS 13.5012 Quite good 
CB-MRAS 7.1828 Good 
GA-EKF 7.2734 Good 
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Figure 12. The change (a) Rotor resistance, (b) The RF-MRAS, (c) The CB-MRAS, and (d) The EKF 
estimated speed 


7. CONCLUSION 

This paper presents three methods of speed estimation of IM: the RF-MRAS, the CB-MRAS, and the 
GA-EKF. Particularly, it focuses on the speedy estimation method using the GA-EKF algorithm with the system 
noise matrix Q and the measured noise matrix R optimized by the GA algorithm. Simulation results in cases of 
variable rotor resistance (constant, decrease, or increase). In the case of constant rotor resistance, the CB-MRAS 
method is much better than RF-MRAS, but the GA-EKF method is also much better than the CB-MRAS. In the 
case of reduced rotor resistance, the CB-MRAS speed estimation method is also much better than RF-MRAS, 
but the GA-EKF method is better than CB-MRAS. With the increased rotor resistance, the CB_MRAS method 
is much better than the RF-MRAS, and the GA-EKF method approximates the CB-MRAS. Thus, in most cases: 
variable or unchanged rotor resistance, the GA-EKF speed estimation method optimized by the GA algorithm 
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which is presented in this paper can claim to be the best, although the CB-MRAS speedy estimation method was 
considered the best in a long time before. 

It can see that the proposed GA-EKF model is less dependent on the induction motor parameters 
than the RF-MRAS and the CB-MRAS model, but the GA-EKF model is the most complex of the mentioned 
models. Therefore, a DSP processor must be used when implementing the GA-EKF algorithm. Today with the 
strong development of semiconductor technology, this method is easily implemented in industrial devices. 
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